Thermodynamics of a Fermi liquid beyond the low-energy limit 
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We consider the non-analytic temperature dependences of the specific heat coefficient, C'(T)/T, 
and spin susceptibility, Xs (T) , of 2D interacting fermions beyond the weak-coupling limit. We 
demonstrate within the Luttinger-Ward formalism that the leading temperature dependences of 
C(T)/T and Xs(T) are linear in T, and are described by the Fermi liquid theory. We show that 
these temperature dependences are universally determined by the states near the Fermi level and, for 
a generic interaction, are expressed via the spin and charge components of the exact backscattering 
amplitude of quasi-particles. We compare our theory to recent experiments on monolayers of He 3 . 
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The Landau Fermi Liquid (FL) theory states that the 
low-energy properties of an interacting fermionic system 
are determined by the states in the vicinity of the Fermi 
surface, and are similar to that of weakly interacting 
quasi-particles. At the lowest temperatures, when de- 
cay of quasi-particles can be neglected, the specific heat 
C(T) oc T and spin susceptibility Xs{T) =const of a FL 
differ from the corresponding quantities for the Fermi gas 
only via the renormalizations of the effective mass and g— 
factor lj. However, this low-temperature limit of the FL 
theory, considered by Landau, cannot tell whether the 
subleading terms in T are analytic, and whether they 
come only from low-energy states (and are therefore de- 
scribed by the FL theory) or from the states far away 
from the Fermi surface. 

For non-interacting fermions, the subleading terms in 
C(T)/T and Xs(T) scale as T 2 and come from high- 
energy states. However, it was found back in the 60s that 
in 3D systems, the leading correction to C{T)/T due to 
interaction with either phonons Q or paramagnons Q is 
non-analytic in T ( T 2 In T) and comes from the states in 
the immediate vicinity of the Fermi surface. The same 
result was later shown to hold for the electron-electron 
interaction kl. More recently, it was shown by various 
groups [ElEullSlilEi that the temperature dependence 
of C(T)/T is also non-analytic in 2D and starts with a 
linear-in-T term. The same behavior was also found for 
the uniform spin susceptibility H, H, S E3 • Two of us 
have argued |(j that, to second order in short-range inter- 
action, these linear-in-T terms originate exclusively from 
scattering of fermions with zero total momentum and ei- 
ther small or near 2kp momentum transfers ( "backscat- 
tering" ) . 

In this paper, we consider the specific heat and spin 
susceptibility for a generic 2D Fermi liquid. In this 
case [if, the leading (constant) terms in C(T)/T and 
Xs{T) are expressed via the two harmonics -Fc^ and 
Fs°^- of the quasiparticle interaction function F(8), or, 
equivalently, via the same harmonics of the scattering 



amplitude A(6) (c and s refer to charge and spin compo- 
nents). We show that the linear-in-T terms in C(T)/T 
and Xs(T) are also universally expressed via the scat- 
tering amplitude, but they are determined by A(9) at a 
particular angle — it, rather than by A(9) averaged over 
the Fermi surface. As there is no simple relation between 
A(n) and F(tt), these subleading terms cannot be simply 
expressed via the Landau function. The only exception is 
a weakly screened Coulomb interaction, when the T-term 
in C(T)/T can be expressed via F(w). 

To shorten the presentation, we discuss in some detail 
the calculation for C(T), and then just present the result 
for Xs(T), which can obtained in a similar manner 
The most straightforward way to obtain C(T) beyond 
the leading term in T is to find the thermodynamic po- 
tential S(T) within the Luttinger-Ward approach [12j, 
and use C(T) = -Td 2 E/d 2 T. The non-analyticity of 
C(T)/T obviously originates from the non-analyticity of 
S(T). The thermodynamic potential 5 is expressed as 
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where So is the thermodynamic potential of the free 
Fermi gas per unit area, Go = (iu>m — efc) , G = 
(ioj m — £fc + £) 1 , £ is the exact (to all orders in the 
interaction) self-energy, and Yj u is the skeleton self- 
energy of order v. Both and the sum S = 
Tj v are evaluated at finite T. Diagrams associated 
with the first two terms in Q correspond to the self- 
energy insertions into the free thermodynamic potential 
S = 2T^ W J d 2 k/ (2ir) 2 lnG (w m , k), diagrammati- 
cally represented by a loop (Fig. ^l a )- O ne can readily 
verify that such diagrams simply renormalize the con- 
stant term in C(T)/T. The non-analytic temperature 
dependence of C (T) /T comes from the third term in . 

To understand the origin of the non-analyticity in 
S(T), consider first a weak short-range interaction /7(g). 
To second order in U, the skeleton term gives rise to di- 
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FIG. 1: Non-trivial second-order and third-order Diagrams 
for the thermodynamic potential. For the Coulomb potential, 
diagrams 1(a), 2(b) and 3(c) represent ring series. 



agrams 2a and 2b in Fig. ^ Assume momentarily that 
U is a constant. Then each of the two diagrams can be 
re-expressed as a product of two particle-hole bubbles 



Tl(q, fi m ), so that 



SE 



cPq 
4tt 2 



n 2 ( g ,n m ), 



(2) 



where 8E = E — So. 

It is intuitively clear that the non-analyticity in SE 
should be related to a non-analyticity in II(<7, fi m ). 
There are two regions of q where II is non-analytic. 
First region is near q = 0, where H(q, D, m ) = 
-(m/27r)(l - \n m \/(n 2 m + {vpq) 2 ) 1 / 2 ). For v F q » 
|O m |, the Landau-damping term (|f2 m |/g) is non-analytic 
in q. This non-analyticity leads to a long-range 
tail of II(r, f2 m ) in real space: H(r,£l m ) oc |O m |/r. 
The second region is near 2kp, where Tl(q,Q ri 
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2kp)/2kF and fl = Q/2kFVF- The singularity at q = 
is known as the Kohn anomaly. Most of the effects 
caused by the Kohn anomaly - Friedel oscillations, Kohn- 
Luttinger pairing - are associated with the (asymmet- 
ric) singularity in the static bubble: Tl(q, ft m — 0) — 
11(0, 0) cx 9 (q) yjq. However, the static Kohn anomaly 
does not cause a non-analytic T-dependence of SE(T) 
in J2J as the momentum integral over static H(q, 0) 
is infrared convergent. The non-analyticity in 5E(T) 
comes from the dynamic Kohn anomaly which is a term 
M/VW\ in n(g,ft m = 0) for -q > \il m \/v F . This 
term leads to a dynamic Friedel oscillation: n(r,f2 m ) oc 
|O m | cos(2k F r)/^/r. 

Integrating in Eq. (0) over the two momentum re- 
gions where II {q, f2 m ) is non-analytic, we find that each 
of these two regions contributes a logarithmic singular- 
ity of the form f2^1n|fi m |, the prefactors are the same. 
This logarithmic singularity is the key effect. Had it 
been absent, the Matsubara sum of fl 2 + ttt, + ... 



would have been controlled by high frequencies, of or- 
der Ep, and would have led to an analytic expansion 
SE =const+T 2 + T A + .... The presence of the logarithm 
changes the story, as now the frequency sum contains a 
universal contribution from frequencies of order T. Using 
the Euler-Maclaurin summation formula, we obtain 
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where u = mU/2ir, and M ~ Ef/2ttT ^> 1. Those terms 
in © that depend on M yield a regular expansion for SB 
in powers of T 2 , whereas the M-independent term - the 
third term in the second line of Eq.© - gives rise to 
SE oc T 3 , and hence to a T 2 — term in C(T). 

Next, we take a more careful look at which four- 
fermion vertices actually contribute to the non-analytic 
part of C(T). For the 2k p -part, the answer follows im- 
mediately from the observation that, for a given direc- 
tion of q, the Kohn anomaly (both static and dynamic) 
comes from the internal fcrmionic momenta near q/2 and 
— q/2 in both bubbles in diagrams 2a and 2b of Fig^ 
The relevant vertex then has the momentum structure 
(k, — k; — k, k), which corresponds to backscattering. 

For the q — part, the momentum structure is less 
obvious as, at the first glance, the internal momenta in 
the two bubbles in diagrams 2a and 2b are uncorrelated. 
However, the logarithmic singularity of the momentum 
integral comes only from the \Q m \/q— part of each po- 
larization bubble. One can show that this part comes 
from the integration over that special region where the 
internal momenta are nearly orthogonal to q. Since the 
relevant momenta in the two bubbles are almost orthog- 
onal to the same vector (q), they must be either nearly 
parallel or nearly antiparallel to each other. We verified 
that the contribution from the near parallel momenta, 
i.e., from forward scattering, vanishes and the full result 
comes from nearly antiparallel momenta. This implies 
that the q — contribution to SE involves a vertex with 
the momentum structure (k, — k;k, — k). This vertex is 
also a part of the backscattering amplitude. 

We can now extend our second-order analysis to a 
finite-range interaction U(q). That only backscattering 
is relevant means that only U(0) and U(2kp) contribute 
to the T -term in C(T) /T. The contribution of the di- 
agram 2b is proportional to U 2 (0) + U 2 (2kp), whereas 
that of the diagram 2a is proportional to U(Q)U(2kp). 
Collecting the prefactors, we obtain for the non-analytic 
part of the specific heat AC = C(T) - jT 



AC (T) /T=-(u 2 + u 2 2kF ~u u 2kF ) 
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where uq — mJ7(Q)/27r and U2k F = mlf(2kp)/2ir. This 
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agrees with the result obtained in Ref. [6( by expressing 
C{T)/T via the self-energy. 

Consider now what happens when we add higher-order 
terms in U. They lead to two types of corrections: self- 
energy corrections to the fermionic lines in the two bub- 
bles and corrections to the four-fermion vertices. The 
self-energy corrections are of the FL type: they account 
for the appearance of the quasiparticlc Z -factors and 
for the replacement of the bare fermionic mass m by 
m*. Vertex corrections generate terms with more bub- 
bles. A generic diagram of n-th order has n bubbles. 
To obtain a T 3 -contribution to 3(T), we need to take 
dynamic, fl m /Q -terms from two bubbles out of n and 
set Q m = 0, Q — ► in the rest n — 2 bubbles, because 
any extra power of 0, m /Q eliminates the logarithmic sin- 
gularity in the frequency integrand in <5S(T). It is in- 
tuitively plausible that once two dynamic bubbles are 
chosen at the n— th order, the rest of the n— th order dia- 
gram constitutes the n— th order correction to the static 
four-point vertex. If this conjecture is true, the series 
of the diagrams for the non-analytic T 3 — term in the 
thermodynamic potential can be re-expressed in terms 
of the two-bubble diagrams in which U(0) and U(2kp) 
are replaced by exact static vertices r(k, — k; k, — k) and 
r(k,-k;-k,k). Accordingly, AC(T)/T is given by the 
same expression as in flljl. but with T instead of U. 

This conjecture, however, needs to be verified as differ- 
ent diagrams for the thermodynamic potentials contain 
different combinatorial factors, and it is a priori unclear 
whether these factors, combined with those counting the 
number of ways two dynamic bubbles can be chosen, 
give the right coefficients in the perturbative series for 
the static vertices. To verify that this is the case, we 
evaluated explicitly the T 3 -term in the thermodynamic 
potential to third order in U(q), and compared the re- 
sult with that given by the two-bubble diagrams with 
the renormalized static vertices, evaluated independently. 
We found that the two expressions are identical. In what 
follows, we assume that this equivalence survives to all 
orders in U(q). 

The renormalization from U(0) and U(2kp) to 
r(k, — k;k, — k) and T(k, — k;— k, k) includes static cor- 
rections coming from the states both away and near 
the Fermi surface (the latter produce powers of static 
II(f2 m = 0, Q — > 0) = — m/2ir U). In other words, 
T(k, — k;k, — k) and T(k, — k;— k, k) include all vertex 
corrections except for the terms coming from the dy- 
namic part of the polarization bubble. In conventional 
notations [l|, t ne fully renormalized T(k, — k; k, — k) and 
r(k, — k; — k, k) are then related to T k (8 — it), which is 
the limit of f2 = and q — > of T(k, p, k + q, p q), 
where both k and p are on the Fermi surface, and 9 is 
the angle between these two vectors. 

The matrix (in the spin space) r fe (7r) is related to the 
quasiparticle scattering amplitude via A(ir) = Z 2 T k (ir) 
[l|. Decomposing the scattering amplitude for a spin- 



invariant interaction into charge and spin components 
A r and A, as 



A(tt) = -^ A c (tt)I + A s {n)a ■ a 



(5) 



and comparing JSJ with a similar decomposition for f fc (7r) 

f fe (7r) = r fe (k,-k;k,-k)J 

-(l/2)r fe (k,-k;-k,k)(i + a-a), (6) 

we obtain 
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^r fc (k,-k,k,-k) = [j4 c (7r) — A s (tt)] ; 

ft]? 



^r K (k,-k;-k,k) = -2— ±A s (n) 

kp 



(7) 



Substituting Z 2 r fc (k, -k, k, -k) instead of [7(0) and 
Z 2 r fc (k, -k;-k,k) instead of U(2k F ) into Q, we find 
the non-analytic part of C (T) in a generic Fermi liquid 

AO(T)/T=-(ffl) (^) 2 [Al(n)+3Al(n)]T. 

(8) 

A similar generalization of the second-order result for 
XAT) i yields 



AXs {T) ^^ F Xs ( ° M ' WT ' 



(9) 



where \s (0) is the spin susceptibility at T = 0. 

Eqs. JSJ) and © are the two main results of the pa- 
per. We see that the non-analytic parts in C(T) and 
Xs(T) are parameterized by two scattering amplitudes, 
A c (tt) and A s (tt). These amplitudes are the new param- 
eters in an extended version of the FL theory, which in- 
cludes non- analytic terms. Notice that A c (w) and A s (tt) 
cannot be simply expressed in terms of the quasiparticle 
interaction function F(0), which, we remind, is related 

to Z 2 f u {6) = (nv F /k F ) (f c {0)I + F s (Q)a ■ a) . The har- 



monics A^ and F^ (a — c, s) are indeed simply re- 
lated, A { a ] = F a (n) /(1 + F { a n) ) in 2D, but A a {ir) is only 
expressed via infinite series of harmonics of the Landau 
function: 



A a W = ^(-irFi"V(l + f| n) ). (10) 

n=0 

Rather involved explicit relations between A a (Tr) and 
F a (n) to third order in U are presented in (l3|. 

Wc find, however, that a simple relation between A(tt) 
and F(tt) does exist if the interaction U(q) is strongly 
peaked at q — 0. In this situation, only corrections due 
to U(0) matter. These corrections come from ring dia- 
grams and can be summed up exactly [the combinatorial 
factor for a ring diagram of order n is (— l) n 2 n ~ 2 (n— 1)]. 
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Evaluating the sum over n, we obtain A c (ir) = F c (tt)/(1+ 
F c (tt)). If, in addition, F s (8) <§; 1, the contribution to 
the specific heat from F c (tt) dominates, and the singular 
term in the specific heat becomes 



AC(T) jT 



3m((3) T 
4tt 



1 + F c (tt) 



(11) 



If F c (ir) is dominated by the n = harmonic, Eq. 
coincides with the result of Ref. 14]. 

The limit F c (ir) — > oo, F s (6) <C 1 corresponds to the 
Coulomb interaction for small r s . In this limit, F c (tt) 
is canceled out from Ijlll) . and the singular term in the 
specific heat becomes 



AC (T) JT 



3mC(3) 
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In agreement with Refs. 0, El> we fi n< l that for the 
Coulomb interaction, the T 2 term in the specific heat 
is independent of r s for small r s . 

Scattering amplitudes in Eqs. JHJ and © can be ex- 
tracted from a measurement of C(T)/T and Xs (T) on 
the same system. To the best of our knowledge, a linear- 
in-T dependence of Xs has not been measured yet. How- 
ever, the linear temperature dependence of C(T)/T has 
been observed in several experiments on fluid monolay- 
ers He 3 adsorbed on graphite 0, 0, 0. To a rea- 
sonable accuracy, the data can be fitted into a form 
C/{NTJE* F ) = -y(T/E* F ), where N is the density per 
unit area in a fluid monolayer, E* F = Epirn/m*) and 
7(2;) a — bx for small x [17| . Both a and b vary some- 
what with N, but the variation is not dramatic, and to 
reasonable accuracy a ~ 3 — 3.3, and 6 ~ 10 — 14 [lfSj . 
According to Eqs.© and 1(8}, a = 7r 2 /3 w 3.3 , and 
& = 0.9[A^ (tt) + 3A 2 (tt)]. A fit to the data then yields 
A\ (tt) + 3^2 (tt) sw 11-15.5 49]. To estimate A c (tt) and 
A s (7r) separately, we assume that the scenario of "almost 
localized fermions" > which describes successfully the 
properties of bulk He 3 , is applicable to the 2D case as 
well. In this scenario, the interaction in the charge chan- 
nel is strong, whereas that in the spin channel is mod- 
erate. A strong interaction in the charge channel means 

(n) 

that Fc 1, in which case the consecutive terms in 

series for A c (tt) [Ea.Q10|)] almost cancel each other, and 
the result is likely to be small. A precsie value for A c (tt) 
depends on how F c decrease with n. However, in two 
model cases Fc = g/(l + n 2 ) and F^ = ge~ n , we 
obtain almost identical results: A 2 (tt) w 0.25 in the 
limit of g ^> 1. This suggests that the observed value 
A 2 (tt) + 3A 2 S (tt) w 11 — 15.5 is almost entirely due to the 
spin part of the amplitude. Neglecting A 2 (tt) , we obtain 
\A S (tt) I « 1.9- 2.3. If the n = harmonic of F s domi- 



nates the result for A 



(tt), i.e.,A s (Tr)^F^/(l + F^ r 
0.7), which is consistent with the 



then F. 



(0) 



3D value F, 



(o) 
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Notice also that if A c (tt) can be neglected compared 
to A s (tt), AC'(T)/T and 5x s (T) contain only one un- 
known parameter (A s (tt)). In this situation, the ratio 
K = AC(T)/(TAxs(T)) is expressed only via the pa- 
rameters describing the leading, analytic parts of C(T) 
and Xs : K = -18((3)m* /tt Xs (0). 

In addition, A s (tt) determines the slope of the linear- 
in-T correction Act(T) to the conductivity of a dirty 
2D FL in the ballistic regime |2l|, which allows to ex- 
press Sxs(T) in terms of Acr(T)"as A Xs (T)E* F /% S (0)T = 
[1 - Aa(T)E* F /a(Q)T] 2 J12. These predictions are 
amenable to a direct experimental verification. 

To summarize, in this paper we showed that the 2D 
Fermi-liquid theory describes not only the leading, con- 
stant terms in the specific heat coefficient C(T)/T and 
the spin susceptibility Xs(T), but also subleading, linear- 
in-T terms. We argued that these terms come from 
backscattering, and are universally expressed via the spin 
and charge components of the scattering amplitude at 
the angle 8 — tt. We extracted the spin component of 
the scattering amplitude from the experimental data on 
C (T) for a monolayer of He 3 . 
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